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Abstract. We discuss the interplay between the piece-line regular and vertex- angle singular 
boundary effects, related to integrability and chaotic features in rational polygonal billiards. 
The approach to controversial issue of regular and irregular motion in polygons is taken 
within the alternative deterministic and stochastic frameworks. The analysis is developed 
in terms of the billiard-wall collision distribution and the particle survival probability, sim- 
ulated in closed and weakly open polygons, respectively. In the multi-vertex polygons, the 
late-time wall-collision events result in the circular-like regular periodic trajectories (sliding 
orbits), which, in the open billiard case are likely transformed into the surviving collective 
excitations (vortices). Having no topological analogy with the regular orbits in the geomet- 
rically corresponding circular billiard, sliding orbits and vortices are well distinguished in 
the weakly open polygons via the universal and non-universal relaxation dynamics. 
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1 Introduction 

Classical polygonal billiards, formed by the piecewise-linear billiard boundary with the vertex 
angles that are rational multiplies of vr, are rational polygons, known as non-chaotic systems 
(see P for review). They are therefore quite distinct from the Sinai billiard (SB) and the 
Bunimovich billiard j3| (BB), in which classical chaotic motion regimes are due to, respectively, 
dispersive effects caused by the boundary, formed by the disk and the square, and the interplay 
between boundary segments, formed by the semi-circles and the square. Moreover, rational 
polygons of m equal sides and m equal vertices (hereafter, m-gons PP) exhibit chaotic-like 
changes in the associated quantum-level spectra |5, but the fluctuations found are very close 
to those known in the universal Gaussian statistics. These controversial evidences for the chaotic 
and regular behavior of polygons are a challenge to dynamic theory of billiards. In the present 
paper we develop a physical insight into the problem, based on analysis of the simulation data 
on the billiard collision statistics. Within a more general context, the delicate problem of the 
interplay between regular and irregular segments, which constitute the total billiard boundary, is 
ultimately related to apparent controversy between the causality and randomness. Our approach 
to the problem is developed through the alternative deterministic and stochastic frameworks, 
which elucidate the duality of chaotic and non-chaotic features in the intrinsic dynamics of the 
rational polygons. 

The singular effects, provoked by the vertex-angle splitting of orbits, is one of the major 
problem in dynamics of non-dispersed billiards pp. Being countable, they cannot be described 
within the standard Markov partition scheme. Similarly to the recent study on the arc-touching 
effects in dispersing SB [5| , we give analysis of vertex-splitting effects on the basis of the genera- 
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lized diffusion equation. The issue of our investigation is the late-time dynamics memory effects 
induced by sides and vertices in both closed and the open rational polygons. 

2 Collision statistics 

2.1 Collision distribution function 

In billiards formed by classical particles of unit mass, moving with unit velocities, the billiard 
collision distribution function D{n, t) is determined as the probability for a given particle to 
undergo n random collisions with the boundary in a time interval t (for the rigorous definition 
of D{n, t) and for the way of its numerical observation, see 0). With accounting for the property 
of ergodicity, this provides the mean collision number 

nc{t) = {n)c= / nD{n,t)dn= / n{x,t) dfj,{x) = — , t » Tc. (1) 

Jo J Tc 

Here n{x, t) is the number of billiard collisions during time t, for a particle of the position set 
X = (x,y,0), which includes its location {x,y) and the velocity launching angle 6 = [0, 27r] 
counted of the x-axis of the billiard table in the space Q. The continues-time evolution of the 
particle orbits preserves the Liouville measure [H] 

dfi{x) = -^^—^dxdydO. (2) 

and establishes the mean collision time (30121 

-Tc = (3) 

controlled by the accessible area A and the perimeter P of a billiard table. Such a description 
of the continuos dynamic system can be derived from the late time behavior of the associated 
collision subsystem |3 |H1 EI] . 

The case of m-gon is specified by circumscribing of the rational polygon below a circle of 
radius r. For a given m, the billiard area Am = (mr^/2) sin(27r/m) and the perimeter Pm = 
2mrsin(7r/m) provide the mean collision time 

TTr / TT\ 

Tern = ^ cos — . (4) 

In the 'infinite-vertex' m-gon limit, one naturally arrives at the circle-billiard (CB) of radius r, 
with the mean collision time 

rcoo = = f . (5) 

A general question arises how geometrical changes of the billiard boundary are followed by 
corresponding classical dynamics. This question is related to the problem of the existence of 
the correspondence principle in polygonal billiards In view of equation ©, it sounds as 

whether the oo-gon is a dynamic analog of the integrable CB. 

2.2 Dynamic diffusion regimes 

We have communicated that the chaotic disk-dispersing SBs and non-chaotic vertex-splitting 
polygonal billiards display a similar behavior attributed to the enhanced diffusion dynamics. In 
polygons, this late-time dynamics is associated with the singular arc-splitting effects, produced 
by singularities of orbits provoked by the vertices pQ. The stochastic approach proposed for 
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Figure 1. Observation of the diffusion dynamics 
in the Sinai (above) and polygonal (below) bil- 
liards. Points: closed circles — numerical data 
for the SBs of side L = 1 for distinct dispersing 
disks of radii R. The open triangle, square, pen- 
tagon, etc. are the data for the regular polygons 
with m = 3,4,5, etc. equal sides. The lines are 
a guide to the eye. The experimental error does 
not exceed the size of symbols. 
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Collision numbers, nln^ 

Figure 2. Observation of the collision distribu- 
tion function. The solid lines are theoretical pre- 
dictions for the CB (jHJ and the 15-gonal billiard, 
estimated in the deterministic-motion approxi- 
mation. The open circles and the point-dotted 
line schematically represent numerical simulation 
of the late-time distributions. For the original 
data see Fig. 1 in |12j . 



analysis of these effects is based on exploration of the generalized diffusion equation, which 
manifests the dynamic relaxation of non-Mar kovian time evolution [TJ. This fundamental 
model-independent equation, which generalizes the Brownian diffusion, is given for the variance 
of particle displacements A^r(t) realized during time t. In application to m-gons, this equation 
was introduced through the variance of the number of random collisions ; namely 



2/z„ 



(6) 



represented with the help of equations and ((IJ). Here stands for the boundary-dependent 
diffusion exponent (which for the Brownian diffusion is equal to 2). In the case of chaotic SB, 
of side L and of disk of radius R, the random walk treatment of the problem provided the 
equation ^ 



A2n(SB)(t)^((n-n(SB))2) 



Tc{R) 



2/z{R) 



t > rc(i?), 



(7) 



corresponding to equation ©, taken with the mean collision time Tc{R), defined in equation 
and the diffusion exponent z{R). Equations © and ((Tj) form a basis for the observation of 
distinct stationary regimes, driven by the billiard geometry. The uniform initial conditions were 
simulated by the standard random number generator RANl^ and the results of the analysis of 
temporal behavior of the directly observed variances for the billiard-wall number collisions are 
presented in Fig. 1. 

^For details of the numerical experiments, see |10|. 
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In SB, the universal (i?-independent) late-time diffusive regime is due to the stabilization of 
particle-displacement statistics with z^^^^^^R) = 1.50 it 0.05. This was rationalized [S] through 
the universal behavior of the particle propagator function ^3] and the Levy flights, which occur 
between the two distant scatterers, belonging to infinite Bleher's corridors of the dynamically 
equivalent Lorentz Gas (LG) [T^. This almost-free motion along the principal diagonal and 
non-diagonal corridors, induced by arc-touching dispersive events, terminates at R > L\/2/4, 
when the diagonal corridor becomes closed. 

In m-gons, the non-universal dynamic regimes are well distinguished of those, observed in 
the chaotic SB and the non-chaotic CB, having the unique ballistic regime defined by z(S^)(0) = 
^(CB) _ -|^^ ^j^g unilateral-triangle and square polygonal billiards, the vertex-splitting effects 
are fairly weak, because the regular orbits are due to the particles mainly move along the infinite 
trajectories, well known in the case m = 3 ^Sl- As seen in Fig. 1, these particles expose their 
ballistic dynamics via z^^^^ ~ z^^^^ ~ 1, characteristic of the integrable CB. In contrast, the 
enhanced diffusion in pentagon and heptagon, developed by the vertex splitting of regular orbits, 
achieves its maximal level, indicated by the highest exponents z^^"^^ ~ z!f^^^ =3/2. In what 
follows, the observed distinct late-time dynamics is examined through the qualitative description 
of long-living orbits. 



2.3 Orbit topology 

It is disturbing that the question of whether every polygon has a periodic orbit remains open 
Meanwhile, any orbit, which is perpendicular to a side returns to this side parallel to itself and 
therefore it is periodic jl()(ll7j. The existence of the perpendicular periodic trajectories in rational 
polygons was proved independently in JH] and ^Hl- The stability of the corresponding periodic 
orbits in arbitrary polygon was also shown |17j . 

In the integrable CB and in 'almost integrable' m-gons, limited by, say, m = 3, 4 and 6, 
the launching angle ^ is almost good integral of motion for the particles randomly injected 
on the billiarcff. Consequently, the orbit-set classification of the late-time living orbits in such 
non-chaotic billiards can be introduced through the collision angle (p = [0,7r/2], estimated from 
the normal to the boundary wall |lfl| I20j . This motion integral provides the basis for our 
deterministic approach to billiard dynamics proposed in J2] and then developed for the square 
billiard [HI and polygons PT| . 

In CB, all regular orbits are ranged between the diameter-beating set and the 'whispering- 
gallery^ sets established, respectively, by the collision angles ip = and 99 ~ 7r/2 (shown in 
Fig. 1). In triangle billiards, the existence of the 'perpendicular orbiV with = was suggested 
in jl6j . and then proved for any m-gon |18l I19j . In the case of square billiard, these orbits 
correspond to the well-known ' bouncing-balV orbits jlOl 1221 05] , given by ~ 0- angle set and 
the if ~ 7r/2-set, moving in the horizontal and vertical directions, associated with the horizontal 
and vertical free-motion corridors. In the even-m-gons, the marginal vertex-to-vertex particle 
trajectories ^ formally correspond to the stable 'diameter-beating' orbits in CB. In contrast, 
the stable sliding orbits, observed for m > 4, are treated as an analog of the whispering-gallery 
orbits in CB responsible for the late-time superdiffusive dynamics, revealed in Fig. 1 in the 
chaotic-like pentagon (2:5'^^'^^ = 1.42) and in the ordered-like hexagon {z^^'^^ = 1.18). The 
chaotic effects are maximal in heptagon {z^j^^^ = 1.48), but with further increase the number of 
vertices, when m > 8, they attenuate gradually. The analysis shown in Fig. 1 suggests that in 
the 'infinite- vertex' limit, the superdiffusive motion is not transformed into the ballistic regime, 
i.e.. Zoo > 1- Although the geometrical correspondence exists, as it is exposed in equation ©, 
the sliding orbits are topologically distinct from the whispering-gallery orbits in the CB, and 
therefore the late-time dynamic correspondence between the 00-gon and the CB does not exist. 
This statement is additionally tested in the next section. 
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2.4 Observation of orbits 

In CB of radius r, the collision distribution function 



°' '<»•') = TS? i-I# ■ "^™' = = (B) 



vrr 



defined in equation is found through the Liouville measure. Similarly, if one ignores the 
vertex-splitting (irregular) effects a given m-gon, the particle-collision function, corresponding 
to the regular orbits, can be obtained through the model (equivalent-side) approximation fH] . 
namely 



Sm (fern 



1 



71 sm ifcm ' ^ 



''"cm 



ii^cm f^cm 

for 99cm cot (^cm < n/ucm < ^cm/ sinc^cm, Otherwise D^""^^ = 0, (9) 

where ipcm = '7r/2m, for the odd m-gon, and (/?cm = vr/m , for the even m-gon. A comparative 
analysis between the distributions observed in 15-gon and CB is displayed schematically in Fig. 2. 

In Fig. 2, the diameter-beating orbit (99 = 0) and the whispering-gallery orbit {ip ^ 7r/2) are 
well observable in CB. In the case of 15-gon, the vertex-singular effects, manifested by deviations 
from the regular-orbit behavior, seem tend to be adapted through the bouncing-ball-like {ip ~ 0) 
and the sliding orbits {ip ~ '?r/2). We see that both kinds of orbits with 99 w 7r/2 (shown by 
open circles and solid squares) are well distinguished, because, the sliding orbits are topologically 
distinct from their prototype in CB. This finding is supported by the distinct diffusion exponents, 
-^15^^^ = 1.2 and z^*-'^) = 1, analyzed in Fig. 1. 

In a given m-gon, the characteristic collision time 

4m^)(^)=r.^ ""^^-"\^°^^"/":\ (10) 

was also found ^211^5 in the regular-orbit approximation. Here il^m = ^cm and V'm = 0, for, 
respectively, the even m-gon and odd m-gon cases, discussed in equation ©. One can see that 
tcmi'p) iS> Tcm for the set 99 « 7r/2. The longest-living regular (sliding) orbits are therefore 
expected to be observed in the weakly open m-gons. 



3 Surviving dynamics 

In the case of weakly open m-gons, one deals with A^(0) = 10^ uniformly distributed point parti- 
cles, which are allowed to escape through a small opening in the billiard boundarj^. A temporal 
behavior of N{t) non-escaped orbits (particles) is commonly scaled by the mean escape time 

11111011111211 

p 

Te = Tc- with A < P, (11) 

where the mean collision time Tc is given in equation ^ . The late-time asymptote of the billiard 
survival probability 

S{t) = ^ oc for t » Te, (12) 



N{0) V t 

describing the algebraic-behavior through the decay dynamic exponent 5. In non-chaotic square 
billiard, the algebraic decay with the exponents 5 ~ 1 was first reported in A careful study 
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Figure 3. Analysis of the temporal evolution 
of the survived orbits in rational polygons with 
small opening width A = O.OSi?. The observation 
times are reduced through the escape characte- 
ristic time To = 300, chosen common for all cases. 
Points: numerical data for the simulated m-gons 
(squares) and the corresponding CB (circles). 



of decay dynamics in the integrable CB and in the 'almost integrable' 4-gon established ^0] two 
distinct channels of slow relaxation. The first channel is due to the regular-orbit motion with 
the decay exponent 6 = 1, provided by the whispering-gallery orbits, and the secondary channel 
originates from the bouncing-ball orbits, exposed ^01 by (5 < 1. 

In chaotic closed and weakly open classical systems (including Hamiltonian systems), ex- 
emplified by BE j23| I24| I25j . infinite-horizon SB j2()[ I26j . and by corresponding low-density 
LG jl4ll27j . the overall algebraic decay was found numerically through the geometry-dependent 
exponents 6 > 1. The channel of algebraic- type relaxation with 5=1, common for both chaotic 
[201 1241 OS] and non-chaotic jlOl I22j billiards, seems to be generic for all incompletely hyperbolic 
systems with smooth convex boundaries. Its independence with respect to the billiard spatial 
dimension its insensitivity to the details of boundary shape JU], including the location of 
a small opening [23], and to the initial conditions ^0] suggests that the late-time a-relaxation 
in classical systems is the universal primary relaxation. Unlike the primary relaxation, the sec- 
ondary relaxation in chaotic billiards with 5 > 1, is shown to be rather sensitive to the billiard 
geometry JU], as well as to the billiard table dimension d and to the initial conditions jlUj . 
Example is the dynamic-exponent constraint 1 < (3 < d estimated j26j and observed in both 
chaotic BB |24j and SB j2()| I28j. As a result, the dynamic-regime classification can be proposed 
for the open classical billiards in terms of the decay dynamic exponent H12]) . namely 



primary channel: 5 = a = 
secondary channel: 5 = (3 



1, universal regime; 

< 1, non-chaotic regime, 
> 1, chaotic regime. 



(13) 



In Fig. 3 we study the decay dynamics in m-gons with m > 8 through the number of survived 
particles ((T^ Nm{t) = Sm.{t)Nm{0). For the particular case of small opening A = 0.05i?, 
the universal relaxation (jl3j) remains stable until m = 64, but when m > m^^^' = 128 the 
primary-relaxation observation window becomes closed. The sliding regular-orbit dynamics 
likely transforms into a singular- orbit chaotic-like motion, when intensively affected by vertices. 



This is indicated by the dynamic decay exponent P-, 



(cxp) 



1.2 and equation (|13() . Also, the 



experimental observations was elaborated for the large opening A = 0.20R, when p., 



(cxp) 



1.1 



and m, 



(exp) 



32 were derived in this case (see Fig. 3 in |21j). We examine that the upper limit 
for the a-channel-observation window shows its sensitivity to the opening width. This suggests 

that the primary relaxation dominates in m-gons with 3 < m < m^ , where the regular-orbit 
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decay motion is established by the universal decay exponent 6m = ex = 1. The latter is due to 
the ^perpendicular^ and sliding orbits, which exhibit, respectively, the non-chaotic and chaotic 
behavior, which are not distinguished within this channel. In order to prove that the observed 
chaotic-like regime in the open system is really caused by the long-living sliding orbits, we deduce 
the observation conditions for such a decay dynamics. 

Similarly to the boundary case R = corresponding to the universal-to-non-universal 

crossover in relaxation dynamics of the chaotic SB (shown in Fig. 1), the a-to-/3 crossover occurs 
in m-gons with m = and is associated with the transformation of the particle motion, induced 
by regular segments of boundary, into that provoked by its singular parts. Just as the periodic 
motion in the crystalline lattice is observed through the collective excitations (phonons), so 
the exploration by particles of translational periodicity of the billiard table with large m can 
provide the collective-particle motion of characteristic time t'^^^\ip) = Tcoo cos~^ {(p) , permitted 
for any regular orbit, preserving the collision angle (p. This qualitative estimate, following 
from equation taken at m ^ 1, indicates that the idealized sliding orbits, observed now 
at ip ~ tt/2 and m ^ 1, leave the piecewise-linear part of the billiard boundary. The real 

sliding-orbit excitations of finite relaxation times rim*^°\ are thought as marginal regular orbits 
established by the maximal collision angles pm^'^^^ = Pm lfTn|> with rim'^°^ = icm"''('/'m)- In 
contrast to the ideal sliding orbits, we also determine the ideal vortex orbits as sliding along 
the piecewise-linear part of the boundary without reflection. Such irregular 'orbits' are formally 
distinguished by the collision angles ifm""^^^ = $m/2, where are rational vertex angles. In 
view of the m-fold rotational symmetry of the billiard table, the existence of the real singular- 
like vortex excitations (or vortices) is expected from the preservation of the angular momentum 
for a certain set of vertex-angle correlated orbits. We therefore assume that the real long- 
living sliding (regular) orbits are precursors of the vortex-like (singular) orbits, attributed to the 
relaxation times rim ^ = 4m'^(^m,/2). In this way, the a-to-/?-relaxation crossover is attributed 
to the observed regular-to-singular sliding-orbit transformation, that occurs starting from the 
large numbers m^^''^ The observation condition for the crossover between the two distinct real 
dynamic regimes can be approximated by the condition of disappearance of the ideal sliding 
orbits above and that of the ideal vortex orbits below ma. 

A small opening of width A (<ti Pm ) can be arranged at any point of boundary [21]. The 
favorable survival situation in the a-relaxation regime, induced by the regular part of the boun- 
dary of length Lm = Pm/m, must exclude the vertex-angle effects under the constraint m < ma- 
The geometrical condition, at which the ideal vortices effectively escape from the billiard table, 
corresponds to the opening of width A ^ L.^, located at one of the angle vertices. Conversely, in 
the late-time /3-relaxation regime, the regular orbits do not survive at all, if any side is absorbed 
by the opening, i.e. when A ^ Lm, with m > rua- Hence, the a-io-(5 crossover in relaxation is 
qualitatively ensured by the condition A = Lm at m = m^. Taking into consideration that in 
polygons with large number of sides Pm ~ Poo = 27rr, one obtains the desirable criterion of the 
a-io-fj crossover namely 

27rr , , 

ma = (14) 

The experimental data m^^^^ = 128 and 32, discussed above, are in agreement with the estimates 
ma = 126 and 31, following from equation H14|l . This finding justifies our surmise that the 
observed long-living vortices in /3-relaxation are due to the sliding orbits modified by intensive 
vertex splitting. Finally, as follows from equation (|14jl in the limit A — > and as experimentally 
supported in Fig. 1, the observation window for vortices disappears in closed polygons. 
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4 Conclusion 

We have rationalized the late-time relaxation dynamics, observed in rational polygons, in terms 
of the so-called sliding orbits, associated with the periodic trajectories. The dynamic properties 
of the corresponding orbits, determined by the collision-angle set close to it/2, are shown to be 
due to the interplay between the piece-line (regular) part of the billiard boundary and its vertex- 
angle (singular) part. The sliding orbits are beyond the simplified polygonal-orbit classification 
by Gutkin where the regular orbits are presented by the 'infinite-past-to-infinite-future' 
trajectories, and the 'mild discontinuity', caused by vertex-splitting effects, is due to the 'infinite- 
past-to- vertex', 'vertex-to-infinite- future', and 'vertex-to- vertex' marginal trajectories. Also, 
they are in a certain way opposed to the perpendicular orbits, rigorously studied in jl7 | I18[ ITO]. 
The experimental prove for the existence of the long- living periodic sliding trajectories challenges 
to further mathematical research in polygons. 

The sliding orbits are well pronounced in simulated dynamics in both closed and open rational 
polygons. When the number of vertices is small (m < 5), these orbits are ineffective in the 
dynamics driving by the long-living bouncing-ball-like orbits. With increasing of the number of 
sides (and vertices), the sliding orbits become stable and therefore regular, but simultaneously 
expose chaotic-like features, in a way similar to that in the chaotic SB. These features have 
been revealed here through the diffusion exponent and the collision distribution function. The 
chaotic-like sliding orbits, nevertheless, play the role of regular whispering-gallery orbits, known 
in the integrable CB. 

In the open m-gons, the chaotic properties of the long-living sliding orbits remain masked 
and hidden, when observed via the primary relaxation within the domain 8 < m < mo,{A). This 
situation changes above the a-to-P crossover, where new collective excitations, denominated by 
vortices, become to be pronounced in the secondary billiard relaxation. These chaotic excitations 
are associated with the collective vertex-splitting effects, established in multi-vertex m-gons, 
which number of vertices exceeds nia = 27rr/A. 

Regardless the proposed rationalization of the regular-like sliding orbits and the irregular-like 
vortices, we have experimentally revealed that the boundary-collision events, observed in the 
collective dynamics in rational polygons are dual of regular and singular of isolated particle 
collisions with the corresponding parts of the billiard boundary. Being related to the zero- 
measure singularities in phase space, these effects violate the integrability of polygons PP, as 
well as the classical-to-quantum correspondence principle |12j . The latter finding is due to 
the memory of vertex-splitting effects, which do not disappear in the billiard dynamics when 
m — > oo, and thereby forbid the interchange of temporal (t — > oo) and the spatial (m — > oo) 
limits. 
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